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FOREWORD 


This  is  the  thirteenth  report  in  the  Data  Reduction  Program  Documentation 
series.  It  is  dated  according  to  the  date  ot  completion  of  the  documentation, 
impiicatio-  is  made  that  this  program  will  not  subsequently  be  modified,  amended, 
or  superceded-  on  the  contrary,  the  history  of  radar  data  processing  is  one  of  con¬ 
tinuous  evolution  of  techniques,  and  it  is  unrealistic  to  assume  that  steady -state  has 
been  reached. 

The  preparation  of  reports  in  this  series  is  under  the  Editorship  of  Charles  R. 
Berndtson  of  Lincoln,  and  of  D.  Ness  man  and  R.  French  of  Phileo-Ford  Corporation. 
Inquiries,  suggestions,  corrections,  criticisms,  and  requests  for  additional  copies 
should  be  directed  to  C.  R.  Berndti  on. 

The  principal  contributor  to  this  report  was  G.  M.  Sheinfeld  (Philco-Ford). 
Due  to  the  intricate,  evolutionary  manner  in  which  the  programs  came  into  being, 
the  editors  regret  that  it  is  in  general  impossible  to  give  due  credit  to  all  --  mathema¬ 
ticians  or  radar  analysts  or  programmers  --  who  contributed  to  the  definition  and 
writing  of  the  programs. 
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( 'QMMON  symbols  and  abbreviations 

(The  units  given  for  certain  quantities  are  the  units  commonly  used  for  those 
quantities,  unless  otherwise  noted. ) 


ADT 

ALCOR 

ALTAIR 

Alt 

APS 

ARS 

ARTP 

ATC 

Avg 

Az 

c 

CADJ 

C-band 

DQLT 

DCO 

El 

EOF 

GMT 

h 

Hz 


ALCOR  Data  Tape 

ARPA -Lincoln  C-band  Observables  Radar 

ARPA  Long-Range  Tracking  and  Instrumentation 
Altitude  (km) 

Average  Pulse  Shape 
ALTAIR  Recording  System 
ALTAIR  Real  Time  Program 
Angle  Track  Console 
Average,  Averaging 
Azimuth  (deg) 


Radar 


Speed  of  Light 

Adjusted  Calibration  Constant  (db) 

ALCOR  frequency,  5664  MHz  (NB)  and  5667  MHz  (WB) 

Wide  Band  Pulse  Doublet 

Designations  and  Communications  Operator 

Elevation  (deg) 

End  of  File 


Greenwich  Mean  Time 

Hours 

Hertz 


IF 

in 

IRV 

LC 

Isb 

min 

NB 

NRTPOD 


Intermediate  Frequency 
Inches 

Inter -Range  Vector 

Left  Circular  Polarization 
Least  Significant  Bit 

Minutes 

Narrow  Band 

Non-real  Time  Precision  Orbit  Determination  Program 


POD  Project  PRESS  Operation  and  Data  Summary  Report 

Phase  Presented  in  deg 

PRF  Pulse  Repetition  Frequency  (pps) 

PRI  Pulse  Repetition  Interval  (s) 

pps  Pulses  per  second 

pts  Points 


R  Range  (km) 

R  Range  Rate  (km/s) 

rad  Radians 

RC  Right  Circular  Polarization 

RCS  Radar  Cross  Section  (dbsm) 

RF  Radio  Frequency 

RGC  Receiver  Gain  Control 

RTC  Range  Track  Console 


s  Seconds 

SD  Standard  Deviation  of  Wake  Velocity 

w 

SDBLT  Wide  Band  Slaved  Pulse  Doublet 

S/N  Signal-to-noise  Ratio 

T  Time 

TAL  Time  After  Launch  (s) 

TGC  Transmitter  Gain  Control 

Tr  Traverse  Angle  (deg) 

UHF  ALTAIR  Frequency;  415  MHz 


V  Velocity 

V  Doppler  Velocity 
d 

V  Mean  Wake  Velocity 

w 

VHF  ALTAIR  Frequency;  155.  5  MHz 

WB  Wide  Band 

WBS  Wide  Band  Slaved 

WTR  Western  Test  Range 

0  Total  Off-axis  Angle  (deg) 

\  Wavelength 

*  Denotes  Multiplication 
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t-LOW  DIAGRAM  SYMIIOI.S 


PROCESS,  ANNOTATION 


DECISION 


CD 


TERMINATOR 


SUBROUTINE:  where  NAME  Is  the  entry 
call  into  the  subroutine 


CONNECTOR:  where  p  specifies  a  page  in  the 
flow  diagram,  and  L  designates 
a  statement  number  in  the  program 
listing  or  a  reference  point  in  the 
flow  diagram 

CONNECTOR:  where  X  implies  a  continuation 
of  the  diagram  to  the  next  page 


INPUT/OUTPUT  OPERATION 


MAGNETIC  TAPE 


PUNCHED  CARD 


DISK 


■ 


altpod 

l.  PURPOSE  AND  UTILIZATION 

A.  Source  of  Data 
ALTA1R1 

B.  Data  Input 
ALTA1R  catalog  tape 

C.  Description 

ALTPOD  is  designed  to  produce  punched  card  metric  data  on  hard  body 
or  chaff  targets  in  a  format  suitable  for  input  to  NRTPOD.  The  data  are  ultimately  used 
to  obtain  a  target  trajectory,  pierce  and  impact  points,  and  a  (3  profile.  ALTP  D  is 
normally  run  every  0.1  s. 

D.  Output 

1.  A  listing  of  computed  quantities. 

2.  Plots  vs  GMT  total  s  of  R  error,  and  Az  and  El  offsets. 
Superimposed  on  these  plots  are  estimates  of  the  standard  deviation  of  R  error  and 

Az  and  El  offsets  based  on  the  S/N  ratio. 

3.  Cards,  punched  in  a  format  suitable  for  input  to  NRTPOD, 

containing  R,  Az,  and  El  corrected  for  biases,  tropospheric  refraction,  ionospheric 
refraction  (optional),  and  angle  offsets  (optional). 
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\ 


Description 


ALTPOD  produces  punched  metric  data  for  input  to;  ■NRTPO&for  both  exo- 
atmospheric  and  re-entry  trajectory  determination.  ALTPOD  automatically  coirects 
the  data  for  known  biases  and  tropospheric  refraction.  Optional  corrections  for 
off-axis  position  and  ionospheric  refraction  are  available.  Ionospheric  refraction 
corrections  should  be  applied  when  the  target  is  in  or  above  the  ionosphere.  Off-axis 


are  listed  to  help  evaluate results. 

The  computations  performed  by  A|/f^OD  depend  on  the  ARTP  version  used  at 
Kwajalein.#  Since  12  March  1971,  R,  Az,  and  Ebare  corrected  for  known  biases 
prior  to  recording  on  the  catalog  tape.  Before  12  March  1971,  R,  Az,  and  FI  were 
corrected  for  biases  by  the  ALTPOD  program. 

A.  R,  Az,  and  El 

R,  Az,  and  El  are  corrected  for  known  biases  and  tropospheric 


refraction  as  follows: 


Prior  to  12  March  1971 


where 


R  =  [I  (85)  +  I  (86)+  1  (3)]  -  AR  +ARt 


I  (85)  +  I  (86)  is  the  unambiguous  range  found  in  FMRDRM, 
R  R 

Items  85  and  86 

1  (3)  is  range  bias  found  in  Calibration  Record  FMRR11, 

Item  3 

All  is  tropospheric  refraction  correction 

AR.  is  ionospheric  refraction  correction  (optional) 


Found  in  FMHDRD,  Items  10  and  11. 


4$ 


(3)  +  Az  offset  (optional) 


Az  =  I  (13)  or  +  I 
■  t  R  v  (t-25  ms)  C 


where 


Az  =  Az  at  time  t 

I  (13),  .  is  Az  encoder  angle  found  in  FMRDRM,  Item  13 

R  (t-25  ms) 

at  time  t  -25  ms 


1^,  (3)  is  Az  bias  found  in  Calibration  Record  FMAACC,  Item  3 


EIt  =  *R  (t  25  ms)  +  lC  ^  "  AE  +  AE.  +  El  offset  (optional) 

where 

EI^  =  El  at  time  t 

I„  (14).  ,  is  El  encoder  angle  found  in  FMRDRM,  Item  14 

R  (t-25  ms) 

at  time  t  -25  ms 

1^,  (5)  is  El  bias  found  in  FMAACC,  Item  5 

AE  is  tropospheric  refraction  correction 

AE.  is  ionospheric  refraction  correction  (optional) 

2.  On  or  After  12  March  1971 


R  =  [I  (85)  +  I  (86)]  -  AR  +  AR.  (optional) 

K  R  i 

Az  =  I  (13)  +  Az  offset  (optional) 

K 

EI  =  I  (14)  -  AE  +  AE.  (optional)  +  El  offset  (optional) 
R  i 


B. 


Az  Offset 


A  , ,  .  VHF  Tr  error  (V) 

Az  offset  (deg)  = - - — 


VHF  LC  sum  (V)  ATr  slope  cos  EI 


*  XK 


where 


VHF  Tr  error  is  found  in  FMRDRM,  Item  19,  if  ALTAIR  is  in  the 
point  target  tracking  mode  and  in  FMRDRM,  Item  99,  if  ALTAIR 
is  in  the  chaff  tracking  mode. 


# 


The  mode  is  found  in  FMRDRM, 


Item  107. 


VHF  Lu  sum  is  found  in  FMRDRM,  Item  23,  if  ALTAIR  is  in  point 
target  tracking  mode  and  in  FMRDRM,  Item  104,  if  ALTAIR  is  in 
chaff  tracking  mode. 

ATr  slope  is  Kund  in  Calibration  Record  FMRR05,  Item  2,  if 
ALTAIR  is  in  point  target  tracking  mode  and  in  Calibration  Record 
FMR5CH,  Items  2-6,  if  ALTAIR  is  in  chaff  tracking  mode. 

XK  is  a  conversion  factor  from  mrad  to  deg 


C. 


El  Offset 


El  offset  (deg)  = 


VHF  El  error  (V) 
VI  IF  LC  sum  (V) 


1 

*  - 

AE1  slope 


*  XK 


where 


t 


1 


VHF  El  error  is  found  in  FMRDRM,  Item  20,  if  ALTAIR  is  in  point 
target  tracking  mode  and  in  FMRDRM,  Item  100,  if  ALTAIR  is  in 
chaff  tracking  mode. 


D. 


AE1  slope  is  found  in  FMRR05,  Item  4,  if  ALTAIR  is  in  point  target 
tracking  mode  and  in  FMR5CH,  Items  12-16,  if  ALTAIR  is  in  chaff 
tracking  mode. 

Ionospheric  Refraction  Corrections 


1.  AR. 

i 

AR.  =  -1.166  (Ry  -  Ru) 

where 

Ry  is  range  from  VHF  range  tracker 
Rjj  is  range  from  UHF  range  tracker 


2.  AE. 


R 


AE.  =  Cot  El  * 

l 


V 


R 


*  4JL.  z 


where 


B  =  R/r  sin  El 
e 

where 


r  is  radius  of  the  earth 
e 


# 


Called  RIOC  and  ELCIO  in  listing. 
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where 


1  +  Alt  (2  r  +  Alt  ) 
_ p  e _ p_ 

(r  sin  El)2 
e 


Alt  is  altitude  of  peak  electron  density 

P  2 

Z  =  1 +  Bcr/(1  +  B  |A  )  *  e  [  1  +  erf  (A)| 

P  ^ 

2 

a  =  a,  /R  sin  El  [  1  +  B  p  (cos  El)  } 
h  P 

where 

a  is  width  of  the  ionosphere 
h 

(JL  =  (  c  -  1)/B 

P 

A  =  (1  -  HLp)/o- 


where 


„  x  2 
2  C  -p  , 

erf  is  an  error  function:  =.  I  e  au 


=  A 


E.  R  Error 


.  .  VHF  R  error  (V)  .  n  v 
R  error  (m)  -  VHF  LC  sum  (vV  K 


where 

VHF  R  error  is  found  in  FMRDRM,  Item  18,  if  ALTAIR  is  in  the 
point  target  tracking  mode  and  in  FMRDRM,  Item  98,  if  ALTAIR 
is  in  the  chaff  tracking  mode. 

VHF  LG  sum  is  found  in  FMRDRM,  Item  23,  if  ALTAIR  is  in  point 
target  tracking  mode  and  in  FMRDRM,  Item  104,  if  ALTAIR  is  in 
chaff  tracking  mode. 

Bm  is  the  range  channel  slope  (yd/V/V).  For  point  target  track  mode, 

B  is  found  in  Calibration  Record  FMRR06  as  a  function  of  track  reference  and  waveform. 
VHF  waveform  is  determined  by  combining  FMRDRM,  Items  3  and  28.  Track  reference 
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(centroid,  leading  edge,  or  trailing  edge)  is  found  in  FMRDRM,  Item  61.  For  chaff 
target  track  mode,  B  is  found  in  Calibration  Record  FMRCHF  as  a  functior  of 
tracking  gate  width  (found  in  FMRDRM,  Item  107).  K  converts  yd  to  ni. 

F.  S/N 

S/N  (db)  =  TGC  +  2.2  -  10  *  (sensitivity  bit) 

where 

TGC  is  VHF  TGC  attenuation,  found  in  FMRDRM,  Item  26.  The 
sensitivity  bit  is  found  in  FMRDRM,  Item  51. 

G.  Standard  Deviations  of  R  Error  and  Az  and  El  Offsets 

The  standard  deviations  (<r)  of  R  error  and  Az  and  El  offsets  are 

computed: 

, ,  v  40  1 

ffAz  offset  (deg)  TSTir*2 

40  _ 

cr_  .  ,,  (deg)  = - *  Z 

El  oftset 

<j  (m)  =  - 

Rerror  TsTn 

where 

Z  is  a  conversion  factor  from  mrad  to  deg 
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III. 


OPERATION 


A.  Input 

Start  and  stop  times  (GMT) 

Sensor  identification 
Date  of  test 
Sampling  interval  (ms) 

Options  for  adding  off-axis  and  ionospheric  refraction  corrections 


A  sample  input 

is  given  in  Appendix  A. 

CARD  1  [(A3,  IX,  312, 

2  (10X,  312,  1X13),  1X14,  2  (4X11)] 

(Col.) 

1-  3 

RADAR 

Sensor  identification 

5-  6 

IYEAR 

Last  two  digits  of  year  - 

7-  8 

IMON 

Month 

Date  of  test 

9-10 

IDAY 

Day 

21-22 

IHR 

(h) 

23-24 

IMN 

(min) 

-  Start  time  (GMT) 

25-26 

ISC 

(s) 

28-30 

ISTl/ 

(ms) 

41-42 

IEHR 

(h) 

43-44 

IEMN 

(min) 

.  Stop  time  (GMT) 

45-46 

IESC 

(s) 

48-50 

IFTH# 

(ms)  J 

52-55 

,  ITV# 

Sampling  interval  (ms) 

60 

IOFF 

Angle  offset  option:  1  = 

add  offsets;  0  =  do  not  add  offsets 

65 

IONC 

Ionospheric  refraction  correction  option;  1  -  add  offsets; 
0  =  do  not  add  offsets 

u 

Must  be  multiple  of  25  ms. 
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CARD  2 
(Col.) 

1-72  LABEL  72  character  label  for  plots 

B.  Output 

LISTING 

GMT  h,  min,  and  s 

R  and  El  ionospheric  corrections  (optional) 

R,  Az,  and  El 
R  error  (m) 

Az  and  El  offsets  (deg) 

S/N  (db) 

Standard  deviations  of  R  error  and  Az  and  El  offsets 
PLOTS 

R  error  and  standard  deviation  of  R  error  vs  GMT  total  s 
Az  offset  and  standard  deviation  of  Az  olfset  vs  GMT  total  s 
El  offset  and  standard  deviation  of  El  offset  vs  GMT  total  s 

PUNCHED  CARDS 


Sensor  identification 

(A3) 

Year 

(3X12) 

Month 

(12) 

Day 

(12) 

h 

(12) 

min 

(12) 

s 

(12) 

Fraction  of  second 

(15) 

Az 

(F9.3) 

El 

(F12.3) 

R 

(F16.6, 

Sample  outputs  are  given  in  Appendix  B. 
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IV. 


PROGRAM  LIMITATIONS 


ITV  <=  9999 

Do  not  add  angle  offsets  unless  an  examination  of  ALTCEP^  shows  they  are 
well  behaved  and  there  is  a  sufficient  S/N. 

Do  not  add  ionospheric  refraction  corrections  unless  target  is  in  range  track 
at  VHF  and  UHF  (see  ALTCEP  listing). 
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V. 


PROGRAMMING 


A.  ALTPDO  (see  Appendices  C  and  D. ) 

ALTPDO  is  the  control  section  of  ALTPOD.  ALTPDO  reads  the  input 
cards,  and  prints  and  punches  the  data.  ALTPDO  calls  CHE  AD,  BMFND,  REA, 

WHDATE,  and  the  plotting  routines.  ALTPDO  also  calls  WHICHV,  a  360-system  subroutine 
indicating  whether  a  job  is  being  run  under  the  time -sharing  (CMS)  or  Batch  (OS)  systems. 


B.  BMFND  (see  Appendix  E.) 

Subroutine  BMFND  contains  a  table  of  constants  used  in  computing 
range  error.  The  call  statement  is  BMFND. 

STORED  IN  COMMON 

B  Range  channel  slope  (yd/V/V) 

M 

C.  CHEAP  (see  Appendices  F  and  G.) 

CHEAD  is  used  to  process  calibration  and  format  records  which  are 
recorded  before  the  data  records  on  the  ALTA1R  transcription  tape.  CHEAD  lists 
format  and  calibration  records  named  in  the  common  statement.  These  are  stored 
and  unpacked  for  later  use  by  the  main  processing  program.  A  sample  CHEAD  output 
is  given  in  Appendix  H. 

The  minimum  size  needed  for  the  item  array  may  be  calculated  by 
the  following  equation: 

Item  size  =  6  *  (Total  number  of  items  stored)  +  (6  *  130) 

Calling  Sequence:  Call  CHEAD  (*) 

*  =  A  return  point  specified  by  a  statement  number 
in  the  calling  program.  Used  for  aborting  job 
by  main  program  if  wrong  tape  is  mounted. 

CHEAD  calls  the  following  subroutines: 

BREADS  (entries  BREADS  and  BREAD);  HDRR  (entries  HDRR 

3 

and  NAMET);  and  FORM.' 
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D.  WHDATF,  (see  Appendix  ].) 

Subroutine  WHDATE  determines  whether  the  ARTP  version  was  on  or 
after  12  March  1971.  The  call  statement  is  WHDATE  (IDARR  (9),  IDARll  (10),  NEW). 

INPUT 

IDARR  (9)  Day  and  first  two  letters  of  month 

IDARR  (10)  Third  letter  of  month  and  last  two 

numbers  of  year 

OUTPUT 

NEW  Date  of  Tape:  1  =  On  or  after  12  March  1971 

2  =  Prior  to  12  March  1971 

E.  REA  (see  Appendices  K  and  L. ) 

Subroutine  REA  reads  and  computes  data  from  the  catalog  tape. 

REA  calls  GET  and  IGET,  HMS,  HM25,  BZERO,  IONCOR,  and  REFC.  The  call 
statement  is  REA  (DTIME,  ICC,  IONC). 

INPUT 

DTIME  Time  of  record  to  be  processed 

IONC  Ionospheric  refraction  correction  option 

INPUT  AND  OUTPUT 

ICC  No.  of  records  processed 

STORED  IN  COMMON 
GMT  Total  GMT  seconds 

RNG2  R 

R  Not  used 

E  El 

A  Az 

OF  Range  tropospheric  refraction  correction 
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F18 

Range  error  (m) 

F19 

Az  offset  (deg) 

F20 

El  offset  (deg) 

N23 

VHF  LC  Sum  (V) 

SNDB 

S/N  (db) 

DTHETA 

Standard  deviation  of  Az  and  El  offsets 

DRANG 

Standard  deviation  of  Range  error  (m) 

GET  and  IGET 

3 

GET  and  IGET  are  entries  to  subroutine  GETS. 

These  routines  will  locate  any  data  item,  unpack  it,  and  interpret  it 
according  to  the  information  in  the  format  table.  They  will  return  the  item  as  a 
binary  integer  (in  the  case  of  IGET)  or  as  a  floating  point  number  (in  tae  case  of  GET). 

GET  (or  IGET)  requires  three  arguments: 

GET  (Format,  Base,  Item) 

Format  Relevant  format  table  address 

Base  Base  address  of  data  block  desired 

Item  Specific  item  number 

G.  HMS  (see  Appendix  M. ) 

Subroutine  HMS  unpacks  the  time  words  found  in  FMRDRM,  Items 
1  and  2,  and  converts  them  into  h,  min,  s,  and  ms.  The  call  statement  is  HMS  (IW1, 
IW2,  IH,  IM,  IS,  IT). 

INPUT 

IW1  FMRDRM,  Item  1 

i\V2  FMRDRM,  Item  2 
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OUTPUT 


IH  Hours 

IM  Minutes 

IS  Seconds 

IT  Decimal  fractions  of  seconds 

H.  HM25  (see  Appendix  N. ) 

Subroutine  DM25  computes  the  TGC  attenuation  found  in  FMRURM, 
Item  25.  The  call  statement  is  HM25  (LALL,  IHMCH). 

INPUT 

LALL  FMRDRM,  Item  25 

OUTPUT 

IHMCH  TGC  attenuation  (db) 

J.  BZERO  (see  Appendix  O. ) 

Subroutine  BZERO  is  necessary  if  a  floating  point  item  is  scaled  BO 

in  a  format  or  calibration  table  description.  BZERO  is  called  after  the  item  has  been 

3 

extracted  by  subroutine  GETS.  BZERO  normalizes  the  data  item  and  puts  the  decimal 
point  in  its  proper  position.  The  call  statement  is  BZERO  (yy). 

INPUT  AND  OUTPUT 
yy  Item  to  be  processed 

K.  IONCOR  (see  Appendix  P. ) 

Subroutine  IONCOR  computes  elevation  angle  corrections  (deg)  to 
compensate  for  ionospheric  refraction.  The  call  statement  is  IONCOR  (ALTP, 
SIGMAH,  RANGE,  ELEV,  DELTAR,  DELTAE). 
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INPUT 


ALTP 

Alt  of  pet.k  electron  density 

SIGMAH 

Width  (km)  of  the  ionosphere 

RANGE 

R  for  each  requested  time 

ELEV 

El  for  each  requested  time 

DELTAR 

VHF  -  UHF  range  difference, 

found  in  FMRDRM,  Item  75 

OUTPUT 

DELTAE  El  correction 

L.  REFC  (see  Appendix  Q. ) 

The  tropospheric  refraction  correction  subroutine,  REFC,  is  based 

4 

on  troDospheric  refraction  tables  in  PPP-36.  A  modified  version  of  this  subroutine 
is  now  in  use. 

The  call  statement  is  REFC  (E,  R,  DEE,  DRR) 

E  Uncorrected  El  (must  be  between  0°  and  90°) 

R  Uncorrected  R  (ft) 

DEE  El  tropospheric  correction 
DRR  R  tropospheric  correction  (ft) 

The  corrected  values  to  be  computed  after  exiting  from  REFC  are: 

El  =  E  -  DEE 
R  (ft)  =  R  -  DRR 

M.  ALTPLT 

ALTPLT  is  the  plotting  routine. 

N.  REREAD 


REREAD  is  a  plotting  system  subroutine. 
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APPENDIX  A 


ALTPOD  INPUT 


ALT  71 0627034432.  000034437.  000  1  000 


340.  0 


150.  0 


CARD  1 

no  000  000  00C000  OOOOOO  00  0000  00000000000  0  0000000  o  oooooooooooooooo 

,  ,  J  4  ,  II  ,  .(Ml  1213  km:  I7I]IS7HI  K»27  »  8  3031  32  33  34  3S  35  W  M  M  40«l  ««  «  «  <S  «M«  <S  5«  51  S2  5]  54  55  56  57  5i  [9  biUI  6!  BJS4  W  67  Cl  tWD  ,  '  ?  1  ' 

mm  Mill  mi  1 1 1 1 1 1 1  m  1 1 1 1 1 1 1 1  iniiimm  min  mum  mmmimimm 
23222222  222222  22222222222222222222222222222222222222222222222222  2222222222222 
3  33333333  33  3  3333  33  3  3333333333333333333333  33  333333333  33333333333333333 

444444444444  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  <»  4  4  4  4  M4 

555556555555555555555555555555555555555555555555555555555555  5  0  5  5  5  5  5  5  5  5  5  3  2551230 

(i  R  6  C  6  E  G  fG  B  6  G  E  C  6  G  B  6  6  G  6  6  6  G  6  6  6  6  B  6  6  6  6  6  6  G  6  6  G  6  6  6  G  6  G  6  G  6  6  6  6  E  G  6  G  6  1  6  6  6  6  6  6  6  6  6  6  E  c  S  6  G  E  G  E  ■  G  5 

1 ! 7 ,  7777  minimum  mmmmmnmmmmmimmmmmmm 
3  5  8  3  ;,  (|  8  C  .1 8  G  6  0  8  8  8  8  8  0  8  8  9  8  8  8  8  8  8  8  8  8  0  8  8  8  8  8  6  0  8  8  8  8  8  8  8  8  8  8  6  8  8  8  8  8  8  8  G  8  8  8  8  8  8  8  8  8  8  r  5 

OS  !  r  q  3  no  f,  n  Ip  0  0  9  0  r  g  9  9  9  9  9  9  0  9  9  9  9  9  <;  3  9  9  9  9  o  0  9  9  9  9  9  9  9  S  0  9  9  9  9  3  ?  9  3 1)  S  8  0  3  o  9  8  8 

iFM  ! 


ALTPOD  PL0T7  FOR  PA  BY  SHEIHFELD 


CARD  2 

n  n  II 0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  3  0  0  0  0  0  0  0  0  0  0  5  J  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0 

!  M  i  m  1 1  m  1 1  m  m  1 1  m  1 1 1  m  m  1 1  m  m m  i  m  i  m  m  m  m  1  m mi  i  m  i  m 

22222222222  22222222  22  22222222222222222222222222222222222222222222222222222222 

3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  <  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  o  3  j 3 

4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  M  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4 

5  5  5  5  5  c  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  r  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  ,i o 

6  G  6  G  G  5  C  6  5  6  6  E  6  6  6  6  G  6  6  6  6  6  G  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  G  6  6  6  6  6  6  6  6  6  G  6  6  6  6  6  6  G  6  G  6  6  6  6  6  6  6  G  6  6  6  E  6  6  6  G  > 

7  j  7  777  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  1  1  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  1 /  7  7  7  7  7 

6  8  8  8  8  8  8  8  6  8  8  8  0  O  8  8  8  8  C  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  0  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  C  8  8  8  8  8  3  8  8  8 
n  o  c  q  o  c  o  qsqqoos  ' p  0  9  9  5  0  9  9  q  9  S  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  S  9  9  9  9  9  9  9  9  9  S  9  9  9  9  9  9  9  9  S  9  099  9  9 

ttM  risen 
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APPENDIX  C 


A  LTPDO  PROGRAM  LISTING 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


(. 

c 


INPUT  PARAMETERS  (ALL  ARE  INTEGERS  EXCEPT  SENSOR  WHICH  IS  ALPHA) 
(AND  PEAK  AND  WIDTH  WHICH  ARE  REAL  FORMAT  F6.2! 
COLUMNS  PARAMETER 

1-  3  SENSOR 

‘>-10  YEAR,  MONTH,  DAY  OF  FLIGHT 

11-IH  START  HOUR,  MINUTE,  STCCND 


17 

18-20  START  THOUSANDTHS  OF  A  SECOND 

21-26  STOP  HOUR,  MINUTE,  SECOND 

27 

28-30  STOP  THOUSANDTHS  OF  A  SECOND 

32-  )5  INTERVAL  IN  THOUSANDTHS  OF  A  SECOND 

<*0  1  =  ACD  OFFSETS,  0  =  NO  OFFSETS  ADDED 

1=USE  I QNS°HER I C  CORRECTION  ON  THE  RANGE,  C=NOT  USED 
50-bS  PEAK  ALTITUDE  FOR  IONOSPHERIC  CORRECTION  FORMAT  F6.2 

60-65  WIDTH  OF  ICNCSPHERE  FORMAT  F6.2 

0=1  SE  TAPE  DATE,  1=N£W  TAPE,  2=GLL  TAPE 
NOTEO  IF  NO  IONOSPHERIC  CORRECTION  PEAK  AND  WIDTH  MAY  BE  C.O, 

Ir  IONOSPHERIC  CORRECTION  AND  PEAK  AND  WIDTH  NCT  SPECIFIED, 
NOMINAL  VALUES  OF  PE AK=  340.0,  AND  WIDTH  =  150.0  ARE  USED 
CARD  2  IS  LABEL  INFORMATION  (72  ALPHAMERIC  LIMIT  FOR  PLOTS)  CS  ONLY 
THIS  CARD  MUST  BE  HERE  AND  NOT  BLANK  FCR  PLOTS 


REAL  *  8  START, STOP, RNG2.DTIME, DINT, XX, GMT, DXT,G3,R3, A3, E3 
INTEGER*2  ITEM 

COMMON/ BE  AD/LN,  IF  LG,  l  ADO , FMRD I D , F MC AT F , F MCS AD , FMCMDB , FMC T I B , EMC  I Db 

1, FMCTD0,FMRDRU,FMRDRM,FMRDRT,FMGL0T,FMRR05,FMAXSP,FMBI AS.FMK5CH 

2 ,  FMRCHF , FMAACC , FMRR 1 1, NAME! 19) ,NI ( 18) ,  IX (  IB  ) , I TEM( 8000 ) 

COMMON  /OKMDTA/GM  T  »RNG2,R,E,A,DF ,F18, f 19, F2C.N23 
COMMON  /SOATA/SNOB, DTHETA.DRAN j 

common  /iocors/rcio.elcio.alti  ,sigh 

COMMON  /KTAPE/NEW 
COMMON  /TITLE/  I  DARR (  1C) 

COMMON  / P4060 /( EG ( 18) 

DIMENSION  G3(2),R3(2),A3(2),E3(2)  ,N0FF ( 2)  ,NTPE (2  ) 

D I  MENS  ION  GP (  ?n0 ),F18P(7C0),F19P(7  00),F20P( 700 ) , SDR (700),SCA(7  00), 
1  SDC  E ( 700 ) 

DATA  NOFF/'NO  • , 'YES  */ 

DATA  NTPF/’NEW  • , ’OLD  •/ 

DATA  ICMS/’CMS  '/ 

DATA  I BL /  ’  '/ 

INPUT  =  5 
IOUTPT  =  6 
N08  =  l 
NO  7=1 

DO  7744  I  =  1,18 
LEG!  I  )  =  I  BL 
7/44  CONTINUE 

CALL  WhICHVI ID) 

IF( ID.NE. ICMS)  GO  TO  10000 
PRINT  9951 

9951  FORMAT ( '  FILE  «,1=YES  ELSE  0,  FILE  7,1=YES,  ELSE  0',/, 

1  '  LNTER  TAPr  DATE  IN  THE  FORM  OF  DDMMMYY*) 

READ ( 5 , 9961N08, N07, I D A RR ( 9 ) , I DARR (  10 ) 
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996  FORMAT  <21  I ,/,  2*4  | 

INPUT  *  ? 

IOUTPT  =  8 
IOOOO  IPP  =  0 

SNOB  =  0.0 
OTHETA  =  0.0 
ORANG  *  0.0 
RC  10  =>  0.0 
ELCIO  =  0.0 
IPLKFV  =  0 
IPS  =0 
LI  *  0 
I  PLOT  *  0 

IF ( IO.NE. ICMSI  CALL  REREAD ( 99 ) 

CALL  BMFNO 
I0R8  *  0 
ICC  =  0 
CALL  CHfcAO 

READ!  INPUT ,99  1 ) R  AO AH , IYEAR,  I  MON, I  DAY , I  HR , I MN , I  SC , I S  TH , 
l  IEHR, IEMN, irsC, I ETH, I  TV, IOFF, IONC ,ALTP,SIGh,MYTAPE 
991  FORMAT (A3, IX 3A7, 2(312, 1)131,1X14, 2(4X11), 2I4XF 6. 2), 4X11) 

IF( 10. EG. ICMSI  GO  TO  995 
MPLOT  =  0 

R  E  AO (  INPUT, 992, EN0=995ILEG 
99?  FORMAT! 18A4) 

00  993  I  ^  1 , 1 A 
I  F  (l.EG(  I  I  .NE  .  I  PL  I  GO  TO  994 
993  CCNTINUE 
GO  TU  995 
9’4  MPLCT  =  1 
995  I  1  CO  *  ISTH/100 

i m i oo  =  hod  *  ioo 
IREST  *  ISTh  -  I M 1 00 

IF( IREST.Gf. 1. ANO. IKEST.LE.24I  ISTh  *  0 
IF ( IREST. GE. 26. AND. IREST. LE. 491  I  STH  »  25 
I F ( IREST. GE. 51 . ANO. .REST.LE.74) I STH  *  50 
IF!  IREST. GE. 76. ANO. IREST. LE. 99)  ISTh  *  75 
IF ( AL TP.LE.0.0 I AL TP  *  340.0 
IFISIGH.LE.O.OISIGH  »  150.0 
NEW  «  MYTAPE 

IF  (  MYT  APE.  £Q.(  I  CALL  WHOA  ,'E  I  I  OARR  (  9  I  ,  I  C  ARR  (  1 0  )  ,  NE  W  ) 

XX  =  I T  V 

OINT  =  XX  *  1.00-03 

XX  *  ISTH 

START  =  IHR*J600+IMN*60*ISC*XX*1. 00-03 
STOP  =  I  tHR*36O0+ I EMN* 60* I E SC* FL OA T ( I ETH) *.C01 
STOP  =  STOP  ♦  7 . 50-02 
PRINT  304.NTPEINEW I 
304  FORMAT!*  *A4,»TAPE') 

WR  I  TE(  6, 303)  IHO,  I,MN,  I  SC,  I  STH,  IEHR  ,  IEHN,  IESC.I  ETH, 

I  I TV.NOFF ( IOFF* 1 ) ,NOFF (  ONC ♦ 1 1 , AL T P , S I GH , 1 0 ARR ( 9 ) , I  0 ARR ( I  0 ) 

303  FORMAT!*  START  * ' 2 (  I  2 , *  0  *  I , I  2 , * . *  1 3 , / , *  STOP  *  ’ 2 ( I  2 , ’ C * ) , I  2 , * . * , 

1  13,/,’  INTERVAL  OF  *14,’  THOUSANDTHS’,/, 

2  *  OFFSETS  AnoEO  *A4,/,*  IONOSPHERIC  CORRECTION  *A4,/, 

3  ’  PEAK  *  *F  6 . 2 , *  WIDTH  *  *F6.2,/,’  TAPE  DATE  *2A4) 

OTIME  =  START 

ICC  *  I 

22  FORMAT ( ’ I  H  M  SEC  RIOC  ELCIC  ELCV  ’, 

1  *  A2  I M  UNAMB  RAN  *,3X’0EL  RA  OEL  AZ  OEL  EL  ’ 

2  *  SNOB  SO  FL  SO  RANGE  SO  AZ*/) 

100  CCNTINUE 
OXT  =  DTIMF 

101  CCNTINUE 

CALL  REAIOTIME, ICC, IONC) 

IF (OXT. LE. 0.0)  OXT  *  OABS(OTIME) 

LI  -  L I  ♦  1 
G3(L I )  =  GMT 
F  3 ( L I )  =  RNG2 
A  3  (  L  I  )  =  A 
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E  3  I  L  I  !  =  E 

OTIMF  =  OTIME  ♦  2.50-02 
IF(LI.EQ.2I  GO  TO  886 
GO  TO  101 
886  GMT  =  G  3 1  II 
RNG2  =  R  3  (  1 1 
A  *  A 3 ( NEW ) 

E  *  E3INEW) 

L I  -  0 

IF  I  1TV.GE.50)  GO  TO  7300 
L!  =  1 

G31  1 1  *  G3( 21 
R3I 1)  =  R3( 21 
A  3 (  1 )  =  A3( 2) 

E3I  1)  *  E3I  21 

7300  IFIGMT.GT.STCPI  GO  TO  200 
OTIMF  =  DXT 
IAFH  *  GMT/3. 60+03 

I  AM  I N * (  GM  T  -  MHH* 3.60+031/6.00+01 
XX  *  GMT  -  IAHH  *  3.60  +  03-1 AMIN*6. 00+01 
I A  SEC  =  XX 

1AHSEC  =  (XX  -IASEC  +5.00-081*1.00+03 
1AHSEC  *  1AHSEC* 1 .00+02 
IAHSCC  =  1 AHSEO / 1  GO 
IF ( I0FF.NE.1 )  GO  TO  538 
E  *  E  +  F20 
A  ■»  A  +  FIR 
538  CONTINUE 

IFIMOOIIPP. 85I.EQ.0.AN0.Na 8.  EC. II  WRITE(IOOTPT,22l 
OEEP  -  OTHETA/ICCSIE/57. 2957811 

IF ( N08. FQ. 1 1 WR I TE (  IOUTPT  ,21 IAHH, I  AMIN, I  AS  EC , I AHSCC , RC I C , E LC I  0 , E , A , 
1  RNG2.F18.F 1 9 , F 20 , SNOB , OTHE TA ,ORANG , CEE P 
IPP  =  IPP  t  1 

I F ( N0  7.EQ. 1 1  WRITE! 7 , 1 07 1 1  RADAR , 1 0R8 , 1  YE AR , I  MON, 10 AY , 

1  I AHH, I  AMIN, I  AS  EC, I AHSEC , A, E ,RNG2 

2  FORMAT!'  '  I  2,  2(  1  XI  2  I  ,  '  .  •  I  3,  1XF8.  3  ,  1  XF  8. 3 , 1  X2  (  F8. 3, 1  X  )  ,  3XF8 . 3  , 

1  3XF8.2,2X2(F6.3,2X),2XF5.1,3XF5.3i8XF6.2,2XF6.3I 

1071  FORMAT (A3,12,1X3A2,3I2,'.'I5,1XF9.3,F12.3,F16.6) 

IPLCT  »  I  PLOT  ♦  1 
GP (  I  PLOT  I  «  GMT 
F18PI IPLOTI  *  T 18 
F 1 9P (  I  PLOT  1  =  r iq 
F2  P( ( PLOT)  *  T20 

soki i plot i  *  orang 

SO A (  IPLOTI  »  OTHETA 

SGCF(IPLQT)  =  0THETA/(C0S(E/57. 29578) ) 

IF ( IPLOT.LT. 700)  GO  TO  5665 

IF(MPLOT.EQ.I)  CALL  ALTPLT(GP,F18P,F19P,F20P,S0R,S0A,SCCE,IPL0T) 
IPLKEY  =  I 
IPLOT  *  0 
5665  CONTINUE 

IFIGMT.GE.STGP)  GO  TO  200 
IPLKEY  *  0 

OTIME  *  OTIME  +  OINT 
GO  TO  100 

200  IGT  =  IGET  ( FMO ATF , I  ADO , 3 ) 

IGS  *  MOC ( IGT, 16 ) 

I GS= I GS  +  (M00( ( IGT/16), 83*101 
I GM  =  MODI ( IGT/256) , 16) 

I GM  *  I  GM  +  (MODI  (  IGT/8096) , 8  I  *  1 0  » 

I GH=MUD ( ( IGT/65536) , 161 
I GH= I GH+ ( M  C  0 ( (IGT/  108 85 76  I, 8)* 101 
IGT=IGH*3600+ IGM*60+IGS 
XIGT  *  IGT 

IF ( ( X IGT  +  IC. 1 .1  E.GMT)  GO  TO  100 
IFIMPLCT.NE.ll  CALL  EXIT 

IF(  IPLKEY. EC. OICALL  ALTPLT(GP,F18P,F19P,F20P,S0K,SDA,SCCE,IPL0TI 

CALL  EXIT 

ENO 
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APPENDIX  D-4 
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APPENDIX  E 


SUBROUTINE  BMFND  PROGRAM  LISTING 


C 


SUBROUTINE  BMFNC 
AUGUST  LC,  1971 
CCMMCN  /BM/  UMM 


BMS (  1,1) 
BMS(  L,  2  ) 


drt  S(  1,3) 
BMS<  2,  l  ) 
3MS( 2*2) 
BM  S ( 2, 3) 
BMS ( 3,  1  ) 
BM  S ( 3,2) 
BMS ( 3,3) 
RETURN 
END 


3,3) 
-2) .687 
-22.937 
-25.313 
-27. B75 
-29.813 
-32. 5CC 
-29.437 
-31.313 
-34.563 


APPENDIX  F 

SIJBPOUTINE  CHEAP  PROGRAM  LISTING 


SUBROUTINE  CHE  AD ( * ) 

VERSION  04/01/71  R.H.  FRENCH 

THIS  SUBROUTINE  IS  USED  TO  PROCESS  ALTAIR  CATALOG  TAPE  OATA 
RECOROS. 

IN  THIS  VERSION  FORMAT  TABLE  FMRR11  HAS  BEEN  AOOEO.  THE  IDARR 
ARRAY  HAS  BEEN  PUT  IN  COMMON  SO  THAT  THE  RTS  AND  ARTP  ^SION 
OATES  CAN  BE  POINTED  OUT  BY  THE  MAIN  PROGRAM.  THESE  CATEISt  CAN 
ALSO  BE  US EO  TP  TRIGGER  THE  FLOW  DF  THE  MAIN  PROGRAM. 


INTEGEK*2  ITEM 

COMMON/BEAO/Ln!iFlG, lAOO.FMROlO,FMCATF,FMCSAD,FMCMDB,FMCTIb,FMClDB 

1 ,  F MCTOB , FNRORD , FMRORM , FMRDRT , FMGLDT , F  MRR05 . FM AXSP , F  MB  I  AS . F MR5CF 

2. FMRCHF.FMAACC.FMRR ll.NAMEI 19 ) . N  I  ( 1 8 ) » I  X (  IB) * ITEM(BOOO) 

COMMON / T I TL E / I P ARR ( 10) 

01  MENS  ION  IBUF1  (  20 4B1.IBUF 2(2048), MCATF(18),NAMEX(19) 

EQUIVALENCE  ( FMRO I D , MC A  TF (  1 ) ) 

EQUIVALENCE  ( I PUF  1  (  1  ) . I BU ) . I  I BUF 2  I  I ) . I BU2 ) 

DATA  NAMEX/'ROtD*.  'CATF'  ,  1 C  SAO ' . 1 CMDB '  .  •  C T I B • , • C I  OB • , ■ CTDB ■ ,  RCRD 
I, • RDRM* ,  'KORT • , 'CLOT • , 'RRUS' ,  '  AXSP' . ' BI AS  •  ,  'R5CH  ,  RCHF  ,  AACC 


2,  'RRL1  '  ,  'HCRC  / 

GAIA  MAX/B192/. IFL/O/ . I T/L/ 
NREC=0 

iJ  =  I 

LU  20  1*1.18 


MCaTFI I ) =0 

f'.A ME  (  I  )  =NAMEX  (  I  ) 

NI (  I  1=0 
I  X  (  I  )  *0 

NAME!  19)=NAMEX(  19) 

CALL  BREADS (LN. I BUF  1 , I BUF 2  .MAX, IFL , I ND X , L EN , 1 FLG , I  ACC ) 
CALL  BREAD ( LN ) 

NREC  =  NREC ♦  1 

I  F (  IFLG.EQ.21GP  TO  55 

I F ( I FLG. EQ. 3. AMO. I T .EQ. 1 ) GO  TO  30 

I  F (  IFLG.EC.31GP  TO  90 

GO  TO  (21,22), INDX 

I  F (  IBU/256.EQ. 1 )G0  TO  90 

GO  TO  34 

IF (  IBU2/ 256.EQ. I ) GO  TO  90 
CALL  NAMET ( I AOO , NAMTD ) 

IF(NAMEO.EC.NAMEX< 19) 1G0  TO  36 

WRITE(6,R5)NAMPD 

FORMAT ( IX , A4, '  FOUNO' ) 

00  18  1=1.19 

tF (NAMED. EC. NA^EX ( I ) ) GO  TO  35 


CONTINUE 
GO  TO  30 

CALL  FORM ( 1  ADO, ITEM! I T) , i B . NAMEO , N TEM , £ 70 ) 
IF(NAMEO.EQ.NAMcX ( 19) )G0  TO  66 

DO  40  1  =  1,  IB 

IF(NAME(  I  ).EQ.NAMED)GO  TO  60 
CCNT INUE 
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GO  TO  30 

WRI  TE!6,56)NREf 

FORMAT!'  PARITY  ERROR  READING  FORMAT  RECORD', 16) 

GO  TO  30 

MCATF!  I  )  =  <B 

WR I TE ( 6 , 9 ) NAMED ,  TEM 

FORMAT!'  FORMAT  = • , AA , ‘  STORED  IN  CCMMCN  N  TEM  = ' ,  I  A  ) 

IX!  I  )  =  I  T 
N I ( I ) =NTEM 
I T  =  I  T ♦6*NTEM 

GO  TO  30 

CALL  HCRR! ( ADD, I TYP, I DARR > 

WRITE<6,72IITYP 
FORMAT!'  TYPE  ',12) 

IF ( I TYP. EG. I )GP  TO  30 
WR I TE ( 6, 68 ) 

FORMAT!  •  TAPE  MOT  CATALOG  TAPE  JOB  TERMINATED  BY  CHEAD') 

RETURN  1 

WRI TE(6,80)NAM^D,NTEM 

FORMAT!'  NAME  =  ' A A , '  NTEM  ='I5,'  *ERROR*  FORMAT  TABLE  LIMITED  T 
.0  700  ITEMS  OR  FORMAT  TABLE  HAS  0  LENGTH') 

DO  81  1=1,18 

IFfMCATF! I ).EQ.0)WRITE(6»82)NAMEX( I ) 

FORMAT! IX, AA, '  NOT  FOUNO’) 

CONTINUE 
WR  ITE<6, 150) 

FORMAT!'  CHEAD  COMPLETE') 

RE  TURN 
ENC 
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APPENDIX  G 

SUBROUTINE  CHEAD  FLOW  DIAGRAM 
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APPENDIX  H 

SUBROUTINE  CHEAD  OUTPUT 


io=as 

r  ype 

1 

uS  A U 

kUUND 

f  OkMA' 

1  =  C  S  A  J 

ST  URf  u 

IN 

COMMON 

NT  EM  = 

t  TIB 

FUUNU 

1 OR  M  AT  =  CT IB 

SrORfcD 

IN 

C  GMMUii 

NTEM  = 

C  1  DB 

kUUNL) 

kUkMAT=C  IUB 

STOKED 

IN 

COMMON 

NTEM  = 

CT  LB 

FOUND 

k  Uk  MA  T  =  C  T  DB 

STOKtU 

IN 

CGMMUN 

NT  EM  = 

KURD 

kOUNU 

k  Ok  M  A] 

[  =  KDRU 

SI  OREO 

IN 

COMM, ON 

N  T  EM  = 

P  D  R  M 

FOUND 

\  OKMAT =  RDRM 

STUPED 

IN 

COMMON 

N  T  E  M  = 

hDKT 

found 

k  HR MAT  =  KDKT 

STOKED 

IN 

COMMON 

NT  EM  = 

1  MSP 

FOUND 

PAC2 

FOUND 

A  PC  3 

FOUND 

A  PC b 

POUNU 

RAO  5 

found 

PA  8 3 

FOUND 

T  KHO 

FOUND 

T  W  MA 

FOUND 

1RTG 

FOUND 

TRMI 

FOUND 

1  KS  p 

FOUND 

KUI  1) 

TuUND 

k  0  k  MA  T 

=  KUI  0 

STOKED 

IN 

CO  1MUN 

N  TuM  = 

A  ACC 

FUUND 

k  Uk  MAT 

=  AACC 

STORED 

IN 

CO  3MUN 

NT  EM  = 

AMP1 

FUUND 

AMP2 

FoUNu 

AMPS 

k  UUND 

AM  P4 

FOUND 

AMP5 

FOUND 

AMPS 

FOUND 

A  PC  l 

FCUNo 

AS  L  P 

FUUND 

asmp 

FUUND 

ATRK 

FUUND 

AXSP 

FUUNU 

I  0  k  N  AT 

~AXSP 

STORED 

IN 

Common 

NT  EM  = 

tC  AL 

FUUNU 

BIAS 

FUUND 

HJPMAT 

=  61  AS 

SI  OPED 

IN 

CUMMON 

NT  EM  = 

BSMC 

FOUND 

L  AT  k 

r-OUND 

F  U  k  M  AT 

=  CATF 

STORED 

IN 

CUMMUN 

NT  EM  = 

CM  AF 

FuUND 

CMbti 

FOUND 

k  OkMA! 

=  CMOh 

stuped 

IN 

COMMON 

NT  EM  = 

GILT 

FOUND 

FUkMAT 

=  GLIJ  T 

STOkED 

IN 

CUMMON 

ntem= 

NOM1 

FUUND 

NO  M2 

FUUND 

N  1M3 

FOUNU 

NO  M4 

FUUND 

UB  J  T 

FOUND 

PHA1 

FOUND 

PH  A  2 

FOUND 

PHA  3 

FUUND 

PHA4 

FDUND 

6 

PH  A  5 

FOUND 

PH  A  6 

FOUND 

3 

KCHF 

FOUND 

FORMAT 

=  RCHF 

STURI  u 

IN 

CUMMON 

5 

R4CH 

FOUND 

R5CH 

FOUND 

4 

format 

=  R5CH 

STOKED 

IN 

COMMON 

RR04 

FOUND 

18 

RR05 

FUUND 

FORMAT 

=  RRC5 

STUKED 

IN 

CLMMUN 

127 

kR06 

FOUND 

KR07 

FOUND 

10 

KRC8 

FUUNU 

RR09 

FOUND 

ERRO 

FOUND 

SCAN 

FuUND 

X  S  EC 

FUUND 

DRG  1 

FOUND 

DRG2 

FOUNU 

DRG3 

FUUND 

SA  SO 

FUUNU 

kRll  NOT  FUUND 

CHEAD 

complete 

4 

8 


X9 

19 

i  i 

5  i 
7 


NT  t  M 


NTtM 
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APPENDIX  J 

SUBROUTINE  WHDATE  PROGRAM  LISTING 


CHOATE 

START 

0 

USING 

*,  15 

STM 

14, 12, 1 2 (  13) 

LM 

7 , 9 , 0  (  1  ) 

L 

2,0(7) 

L 

3,0(8) 

LR 

4,3 

SRDL 

4,24 

SRL 

5,16 

C 

5,=X'00C0F7F1* 

6L 

OLD 

BH 

NEW 

SRDL 

4,8 

LR 

4,2 

SRDL 

4, 16 

SRL 

5,8 

A 

5, =X' 4 0000000' 

LA 

6, TABLE 

LA 

7,0 

TEST 

C 

5,0( 7,6) 

BE 

LOOK 

A 

7 , =F ' 4  • 

B 

TEST 

LOCK 

C 

7,=F'fl* 

BE 

Til 

BL 

OLD 

NEW 

L 

3 ,  -F  '  1  • 

RET 

ST 

3,0(9) 

LM 

14,12,12(13) 

BR 

14 

OLD 

LA 

3,2 

B 

RET 

T)  I 

SRL 

2,16 

C 

2 , =X ' POCOF  IF  1  ' 

BL 

OLD 

B 

NEW 

CNOP 

0,8 

TABLE 

DC 

CL4  *  JAN* 

DC 

CL4  *  FEB' 

DC 

CL4  *  MAR' 

DC 

CL4 '  APR' 

DC 

CL  4 '  MAY' 

CC 

CL  4 '  JUN' 

DC 

CL4 '  JUL' 

DC 

CL4 '  AUG' 

DC 

CL4 '  SEP' 

DC 

CL4 '  PCT* 

DC 

CL4 '  NOV' 

CC 

LTORG 

END 

CL4 '  PEC' 
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APPENDIX  K 

SUBROUTINE  REA  PROGRAM  LISTING 


SUBROUTINE  R E A  I OT I  ME ,  ICC, I ONC ) 

C  AUGUST  10,  ini 

REAL  «  B  START, STOP, RNG2.0TIHE, DINT, XX, GMT 
IN  TEGER*2  ITEM 

COMMON  /0RM0TA/GMT,RNG2,R,E, A , OF , F 1 8,  E 1 9,  F20, N2 3 
COMMON  /S0ATA/SN08,0THETA,0RANG 
COMMON  / IOCORS/RC 1 0 , EtC IC , Al TP , S IGH 

C  COMMON /BE AC /LN, I FLG , IAOO  ,FMR0 10 .FMCATF , FMCS AC , FMCMDB , FMCT 1 B , FMC I DB 

1 , F  MCI  OB, FMRORO, FMRORM, F MR  CRT .FMGLOT ,FMRR05, FMRR06 , FMRR 1 1 , F MR 5CF 

2, FMAP0l, FMAP03 .FMRCHF , FMAACC.NAME ( 20  I , M ( IP  I , I  X ( IP  I , 1 T EMI BCDO ) 
COMMON/ BE AO/LN, I FLG,  1  ADO , FMRO 1 0 , FMCATF , FMCS AO , FMCMD B , F MCT 1 B , FMC I  OB 

1, FMCI0B,FMR0R0,FMR0RM,FMR0KT,FMGL0T,FMRR05,FMAXSP,FMBIAS,FMR5CF 

2,  FMRCHF, FMAACC, FMRR  11 .NAME!  19) ,NI I  IB) , IX  I  18 1 , 1 TEMl  8000 1 
COMMON  /BM/8MS(3,3) 

COMMON  /KT  APE/ NEW 

DIMENSION  IBUF 1 1  2048) , 18UF2 I 206B > , MC A TF I70),NAMFXI70>, BMC  I  5) 

I NTEGER*2  I BU 

EQUIVALENCE  (FMRO 10, MCA  IF  I  II) , I IBUF 1 1  1 ) ,1  BU) 

OAT  A  NAMEX/'ROIO' , ’CATF’ , 'CSAO1 , ' CMOB *,*CTIB','CIDB,,,CTOB,,,RCRC’ 

1 ,  • ROHM ','RORT', ’FMS,','APOl’, 'PA02* , 'AP03,,'PA05>,,PA6j' 

2 ,  ' TRHO '.'TRMA'.’TRl  • , *  TRHI  •  , * IRSP*/ 

OATA  SLAP/1. 0H1 1/ 

OATA  BMC /-1966. 8,-3893.7,-7787. 6, -155 75. 0,-31169.0/ 

DATA  1TIM/1/ 

ICATIlML)  *  IGF T( FMCAi F, IAOO, IML) 

SORMIIML)  -  GET  I FMRORM,  IBASM, IML I 
IOP.MIIML)  =  IGFT  I  FMRORM,  IBASM,  IML) 

FIB  =  0. 

F 1 9  *  0. 

F  20  *  0. 

N23  *  0 

GO  TO  I  ID, ICO) , I  TIM 
10  CALL  BREACILN) 

ICC  *  ICC  *  I 
IF!  IFLG.EQ-1  1  '‘O  TO  66 
GO  TO  10 
66  ITSI  «  1 

I  TSJ  *  1 

N50  *  ICATH2) 

L  60=  ICATI  13)  *  3 
(BASE  «■  ICATI  191  *  3  ♦  IAOO 
C  I  BASE  »  I  BASE* I  IGET(FMR0R0,I8ASE,l)-l)»3 

100  00  106  I  *  ITSI.N50 

LORM  =  I  BASE  ♦  I  , - l I  *  L 50 
00  103  J  =  I 'SJ , 2 
IBASM  *  LORM  ♦  (J-l)  «  120 
I  W  1  =  I  OR  M  I  1) 

IW2  =  :i)RM(  2) 

CALL  HMSI I  HI , IW2, IH, I M, IS, I T I 

I F I OT IME.LF.0.0 )  OTIME  *  I  IH* 3600* IM*60+ 1 S* F*  CA I  I  1 T ) * . 00 1 ) 

KL  =  OTIME 

XP  =  I  h* 36C0  ♦  I M * 60  ♦  IS 

GMT  =  KP  ♦  FLOAT! IT)  ♦  0.10-02 

IF (KP.LT.KL I  GO  TO  103 

IF  I  GMT .L T .CT IMT )  GO  TO  103 

ITIM  *  2 

R  *  RORMI 1C) 

ACOR  =  0.0 
ECCR  =  0.0 

IF  I  NEW. EQ. 2 )  ACUR  *  GE T I F M A  ACC , I  ADO , 3 )  »  57.29578 
IFINEVi.EC.2)  ECOR  =  GE  T  I  FMA  ACC ,  I  AGO  ,  5 )  *  57.2957B 
A  *  RORMI 13)  *  360.0  ♦  ACUR 
E  =  RORMI 16)  ♦  360.0  ♦  ECOR 
COSEL  *  1.0/CCJr  I  E/57.2957B) 

F9  «  I ORM I  9 ) 

ICHAFF  *  I  OR  M  I  108) 

IF!  ICHAFF.EO.il  GO  TO  21 
F 1 8  *  RORMI 18) 

F  19  =  RORMI  19) 

F  20  =  R0RMI20) 
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N2J  *  10RHI23) 

IFIN23.EC.0I  GC  TO  27 
I T  20  *  IORHI 28  ) 

I  T  6  I  *  I0RM(6LI 

IF(  IT28.GT.A.0R.IT61.GT.AI  GO  TO  20 
IF  (  I  T6  l.fcO.  A  I  I  T6  1  =  3 
I F ( IT28. EO.  A  I  ! T28  =  3 
BTHING.*  BHSI  I T28, 1 F6 I  )  *  .91AA 
F 1 8  *  (F18/N23)  *  BTHING 

20  XA  *  GE  FIFHRR05, | AOC.A I 
X2  *  GEFIFHRR05, [ACC, 21 
CALL  UZEROI XA  ) 

CALL  8ZER0IX2I 

F 1 9  «  (F19/N23I  *  1.057/X2I  *  CCSEL 
F 20  «  (F20/N23I  *  I.057/XA) 

GO  TO  2  7 

21  ITTH  =  IORHI  1071 
USEBH  *  BHCI ITEM) 

F98  =  IORHI 98  I 
F10A  *  IORHI 10A) 

IFIFIOA.EC.O.O)  GO  TO  27 

FI8  *  (F9E/F10A)  *  USEBH  *  .9IAA 
F99  »  IURNI99I 
F I  00  *  IORHI 100 ) 

I  107  *  IORHI  107) 

IA  *  1107  ♦  1 

IE  =  1107  ♦  II 

XA  =  GET  I FHR5CH, )  ADC  t IA) 

XE  =  GET! FHR5CH, IADC, IE ) 

XA  *  BZERCIXA) 

XE  *  B2ER0IXE) 

F 1 9  =  (F99/F10AI  *  I.057/XA1  *  CCSEL 
C  F 1 9  =  IF99/F10A)  *  (0.05  7/0.053)  *  COSEL 

F20  *  (F1C0/F10A)  *  (.057  /  XE ) 

C  F20  *  (HOO/riOA)  *  (.057  /.053) 

27  FB5  ~  IORHI 85  I 
F 86  *  R0RHI86) 

RNG2  =  I F  8  5 ♦ F  86 i  *  3.0 
KCIC  *  0.0 
ELCIO  *  0.0 

IF! IONC.EQ.OI  00  TO  62 

OTR  -  WORM  I  75)  *  3.0 

RCIC  «  OTR  *  1.166  *  .30A8E-03 

OTR  *  OTR  *  . 30A8E-C3 

RKH  *  RNG2  *  .30A80-03 

EEE  *  E 

CALL  IONCOR(ALTP.SIGH,FKH,cEE.DTk,ELCIC) 

6’  COR  »  0.0 

IFINEW.EQ.2)  COR  -  IGET(FHRR11,IAD0,3)  ♦  1.87371 

F  *  RNG2  ♦  COR 

CALL  REFCIE.F.CE.OFI 

RNG2  »  (IRNG2  ♦  COR  -  OF)  *  .30AB0-03)  -  RCIC 

E  »  E  -  CE  -  El  CIO 

SNOB  *  0.0 

OTHETA  *  0.0 

ORANG  =  0.0 

I01CT  =  IORHI 26 ) 

CALL  *  I  CRH  I  25 ) 

CALL  HH25ILALL. IFHCH) 

SN0B*2.2<- ICI0T*3- 10*10RH( 51 )  ♦  IHHCH 
PRF  =  1.0  /  I  1C  AT  I  5 )  *  1.0E-06) 

SN  *  10.0  **  I  SNOB / 10. 0 1 
SS  *  SCR  T I SN  ♦  0.025  *  PRF) 

OTHETA  =  (AO.O/SS)  *  0.05729578 
ORANG*  21.75/SS 
277  ITSJ  =  ITSJ  ♦  1 

IFIJ.GE.2)  GC  TO  101 
GO  TO  102 

101  ITSI  =  ITSI  ♦  I 
ITSJ  *  1 

IF! ITSI.GT.N50I  ITIH  *  1 

102  RETURN 

103  CONTINUE 
10A  CONTINUE 

GO  TO  10 
ENO 
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APPENDIX  L 

SUBROUTINE  REA  FLOW  DIAGRAM 
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APPENDIX  L-2 


DETERM’NE 
ADDRESS  OP 
MRDRM  BLOCK 
2  EE  PROCESSED 


JN  PACK 
"iME  OF 
FMRDRPl 
B  LOCK 


CONVERT 

"I ME  F ROM 
R  M  IN,  S,  t  MS 
TO  GMT 


>  DTINE 


A"Z-  $  EL 
BIASES  “0 


APPENDIX  L-3 


* 


^VEf?£,ON  ON 
\CR  AFTER 

NlA  3  / 
\71/ 


*C 


GET  A 3  £  &L 
31  AGES  FRCMl 
FMAAOC. 
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APPENDIX  L-5 


CCMPLTE  S  N t 
S  TANDARD  DEVI  - 
A',OKi  OF  R  ERROR 
£  A2  5  EL 
OFFSETS 


/fMRDRM  \ 
BLOCKS  IN  RECOR 
XpROCESSEDv'' 


z  FtARDRh 
sB  LOCKS  IN  RECOi 
XpPOCESSED/ 


APPENDIX  M 


SUBROUTINE  HMS  PROGRAM  LISTING 


SUBROUTINE  HMS I  I W 1 , l W2 , I H , I M ,  l  S ,  I  T  ) 
MU  =  MOD ( I Wl« 16) 

MT  =  MGD ( (  I W 1 / 1 6 ) t  8  ) 

IM  *  MT  *  10  ♦  MU 
I HU  =  MOD ( ( IW1/256) , 16) 

IHT  *  MIDI  I  IW1MC96)  ,4) 

IH  =  IHT  *  1C  +  I  HU 
I T  TH  =  MODI IW2« 16) 

ITHU  =  MODI  I  IW?/16),16) 

ITTN  =  MODI  I IWP/256) • 16) 

IT  =  I T  TH  ♦  ITHU  *  10  +  ITTN  *  100 
ISU  =  MODI  I IW2/4C96) , 16) 

1ST  =  MUDI I  IW2/65536) ,8  ) 

IS  =  ISU  +  1ST  *  10 

RETURN 

END 
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APPENDIX  N 


SUBROUTINE  HM25  PROGRAM  LISTING 

SUBROUTINE  HM25 { L ALL » I HMCH ) 
DIMENSION  K { 5 ) 

IHMCH  =  0 
I  ALL  =  LALL 
DO  3  I  =  I, 

K(I)  =  0 

K  (  I )  =  MODI  I  ALL  »  2  J 
I  ALL  =  I  ALL /  1 0 
CONTINUE 

IFIK(l).EC.l)  K(l)  =  3 

IF(K(2).EQ.l)  v{ 2)  =  3 

IF(K(3).EQ.l)  *(3)  =6 

IF(K(4).EQ.l)  =  ip 

I F  ( K  (  5  )  .  EQ.  1)  K  {  5  >  ^  24 

IHMCH  =  K{ l)+K(2)+KI3)+K(4)+K(3) 

END 


APPENDIX  O 

SUBROUTINE  BZERO  PROGRAM  LISTING 


START  0  VERSION:  7/29/70 

ENTRY  BZERO 
USING  *,15 


bZERO 

B 

H5 

CC 

X  '  05  '  , CL5  *  BZERO  • 

H5 

STM 

14,12,12(13) 

L 

6,0(1) 

M  VC 

WCRC,0{ 6 ) 

LF 

4  ,WCRP 

L  PER 

6,4 

CE 

6, =E  *  1  .  • 

BL 

H20 

AU 

4  »  X6 

STE 

4  ,  WCRP 

l 

7,WCftT 

SR 

9,9 

LTR 

7,7 

SLL 

7,8 

SNL 

H  1 0 

LA 

9,1 

H1C 

SRL 

7,8 

SR 

2,2 

AH 

2, =X ' 01C3  * 

SRDL 

2,2 

SRL 

3, 3C 

SLL 

7,0(3) 

SRL 

7,2 

SRDL 

2,8 

OR 

7,3 

SLL 

9,31 

OR 

7,9 

SER 

2,2 

ST 

7 , WCR  0 

AE 

2 , WCRP 

STE 

2,WCRr> 

MVC 

0(4,6) , WORD 

H20 

RETURN 

(14,12) 

WORD 

DS 

IE 

X6 

DC 

X* 46000000' 

END 
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APPENDIX  P 


SUBROUTINE  IONCOR  PROGRAM  LISTING 


SUBROUTINE  I ONCGR ( AL T P , S I GM AH , R ANGE , E L E V , DEL T AR , DEL T AE ) 

C  ALTP^ALT I TUHE  OF  PEAK  ELECTRON  DENSITY  (INPUT  IN  KM) 

C  SIGMA H  =  WICTH  OF  I  HE  IONOSPHERE  (INPUT  IN  KM) 

RANGE  =R ANGE  FOR  EACH  REQUESTED  POINT  FROM  MAIN  PROGRAM  (INPUT  IN  K 
C  ELEV  =  EL£VATT  ON  FOR  EACH  REQUESTED  POINT  FROM  MAIN  PROGRAM 

C  l( INPUT  IN  DPGRRES) 

C  CEL T AR=RDRM ( 75 >  /1093.611 

C  ELCU (^CORRECTED  EVELVAT ION  ANGLE  FOR  EACH  REQUESTED  POINT  (OUTPUT 

C  1  IN  DEGREES) 

ELEV=ELEV/57. 29578 

ERAD  =  6378.  l'«5  VERSION:  5/13/71 

S I NEL=S I N ( ELEV ) 

BETA  =RANGE/(ERAD*SINEL> 

C  =  1+ (  ALTPM  2*ERAD*ALTP)  )/ (  (ERAC*S  INEL  )  **2) 

SQRTC  =  SGRT( 0  ) 

UP=<  5QRTC-1 ) /BETA 

5  I GMA  =  S I GMAH/ ( RANGE*SI NEL* ( 1  +  BE  T A*UP* ( CC S ( E LE V ) * *2 )  ) ) 

A  = (  1-UP  >/SIrMA 

2  =  1+ (  ( BET  A*S IGMA  )/ ( 1  +  BETA*UP ) )* (EXP (- ( A*A ) ) / ( . 8862269* (  1  +  ERF ( A  )  )  )  ) 

CEL TAE= (COT  AN ( ELEV) M  DEL  TAR /RANGE )*<(1+BETA)/C))*Z 

CELTAE=CELTAE*57. 29578 

RETURN 

END 


48 


APPENDIX  Q 

SUBROUTINE  REFC  PROGRAM  LISTING 


100 

120 

200 

220 


300 

320 

34C 

400 


SUBROUTINE  REFC  (  E , 5 ,  DEE  ,  ORR  ) 

01  HENS  ION  DEI  16,8  / ,  OR  I  16,  8  1  ,EDI  16)  ,RDt  8)  VERSION:  6/16/70 

DATA  OE/O.O  ,0.0  ,0.0  ,0.0  ,0.0  ,0.0  ,0.0  ,0.0  , 

10.0  ,0.0  ,0.0  ,0.0  ,0.0  ,0.0  ,0.0  ,0.0  ,0.0313, 
20. 030 3,0. 0292,0. 028 7, 0.02 82, 0.02 72, 0.0262, 0.0253, 0.0243, 0.0223, 
30. 0214, 0.0195, 0.0171, 0.0135, 0.0075, 0.0  ,0.09 37, 0.0848, 0.07 70, 

40. 0732, 0.0694, 0.062 7, 0.05 7 1,0. 0522, 0.04 80, 0.0412, 0.0385, 0.0 337, 
50.0278,0.0205,0.0105,0.0  ,0.1850,0.1520,0.1250,0.1140,0.1050, 

60.0904,0.0795,0.0  708,0.06  36,0.0523,0.0478,0.0405,0.032  3,0.0  2  29  , 
70.0114,0.0  i 0.53 10,0. 3070, 0.2  120, 0.  1 830, 0.1600, 0.  1280,0.  1060, 

80.0899,0.0780,0.0612,0.0550,0.0455,0.0354,0.0246,0.0120,0.0  , 

90.7550,0.3720,0.2400,0.2020,0.1750,0. 1370,0.1120,0.0942,0.0811, 
40.0631,0.0566,0.0466,0.0361 ,0.0250,0.0122,0.0  ,0.9120,0.4110, 

80.2560,0.2140,0.1840,0.1420,0. 1 1 50 , 0 . 096 7 , 0 . 08 30 , 0 . 064 3 , 0. 0575 , 
CO. 0472, 0.0 365, 0.0252, 0.0122, 0.0  ,0.9700,0.4200,0.2600,0.2200, 

00. 1900,0. 1460,0. 1 1 70,0.0980,0.0840 ,0.0653,0.0584,0.0478,0.0369 , 
E0.0254, 0.0123, 0.0  / 

OATA  OR/  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0, 


1  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  C.O,  22.6,  21.5, 

2  19.4,  18.5,  17.6,  16.8,  16.1,  14.8,  14.2,  13.2,  12.0, 

3  7.7,  67.3,  57.9,  50.2,  47.0,  44.1,  39.3,  35.4,  32.1, 

4  22.9,  19.7,  16.3,  12.7,  9.4,  8.1,132.0.  98.5,  77.4, 

5  52.9,  44.7,  38.4,  33.4,  26.4,  23.9,  20.1,  16.4,  12.7, 

6340.0,167.0,103.0,  86.1,  73.4,  56.7,  46.2,  38.9,  33.6, 

7  20.2,  16.4,  12.8,  9.5,  8 . 2 , 405 . 0 , 1 7C . 0 , 1 04 . 0 ,  86.3, 

8  46.3,  38.9,  33.7,  26.5,  24.1,  20.3,  16.5,  12.8,  9.5, 

9171.0,104.0,  86.6,  73.9,  57.1,  46.4,  39.0,  33.8,  26.8, 

A  16.6,  13.0,  9.8,  8.4,446.0,172.0,105.0,  87.4,  74.0, 

8  39.2,  34.0,  27.0,  24.6,  20.7,  16.7,  13.0,  10.0,  8.4/ 

DATA  ED, RTDEG/9. 01, 2.0, 4.0,5. 0,6. 0,8.0, 10.0, 12.0,14.0, 
124., 30., 40., 60., 90., 5 7. 295 73/ 

OATA  RD/0.01, 10. ,30., 60. ,200. ,400. ,1000. ,2000./ 

IFIR.LE.0.01G0  TO  300 

RG-R/6080.27 

00  100  IED-2,14 

I* 17-1  ED 

IFIE.GE.EDI  I  1  ICO  TO  120 
CONTINUE 


20.4,  19.9 

10.4,  8.6 

29.3,  24.8 
69.7,  63.2 

9.4,  8.1 

26.4,  24.0 
73.6,  56.8 

8.2,421.0 
24.3,  20.5 
58.0,  46.6 

18.,— '0., 


1*1 

DD  200  JRU*2 , 8 
J* 10- JRD 

IFIRG.GE.ROI J11G0  TO  220 
CONTINUE 


J*1 

IFIJ.EQ.8IG0  TD  340 

ZR=ALOG( RG/ROI Jl ) /ALOG I RD I J* l ) /rq I J ) I 

IFIE.LE.O.O) 30  TC  320 

ZE=ALUG( E/EDI  I) )/ ALOG I  ED  I  I ♦ 1 ) /ED  I  I  )  ) 

DEI* I (DEI I+1,J)-CE( 1, J))*( 1.-ZRIMDEI 1  ,  J+l  l-DEI 1, J))*ZR)*ZE 
DE2*( (DEI  I , J+11-0EI I, J  »  » *  I l.-ZE)* (OEI 1*1, J+ll-DEIl , J  +  l ) )*2EI*2R 
DEE*DE1+DE2+DE( I , J ) 

OR  1*1 (OR! I *1, Jl-DRI 1, J) )*l 1.-ZR1+ ( DR  I  I  ,  J+l I -OR! 1, J)  )*ZR)*ZE 
DR2* I  I  OR  I  I , J* 1 l-ORI 1,J)>*(1.-ZE)+(DR(I+1,J+1)-DR(I,J+1))*ZE)*ZR 
DRR=I DR1+DR2+DR I  I , J ) ) 

GO  TO  400 
DEE*  0. 0 
0RR*0. 0 
GO  TO  400 

DEE*Ce(I,J)+(DF((,J+l)-DE(!,J)) *ZR 
DRR*DR(I,J)+(DR(l,J+l)-OR(l,JI)*ZR 
GO  TO  400 

DELT-IE-EOI I  1 l/IEDI l+ll-EDI 1)1 
OEE*OELT*(DE( 1+1, JI-DEI 1 , J)  )  +  0E I  I , J) 

ORR*DEL  T* I  DR  I  I ♦ 1 , J 1 -OR  I  I , J ) ) ♦ DR  I  I , J ) 

RETURN 

END 
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